What do we (not) know theoretically about solar neutrino fluxes? 
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Solar model predictions of ®B and p-p neutrinos agree with the experimentaUy-determined fluxes 

(including oscillations): 0(pp)measured = (1.02±0.02±0.01)(^(pp)thoory, and 0(*B)measured = (0.88± 

0.04±0.23)(^(*B)thcory, lo" experimental and theoretical uncertainties, respectively. We use improved 
input data for nuclear fusion reactions, the equation of state, and the chemical composition of the 
Sun. The solar composition is the dominant uncertainty in calculating the **B and CNO neutrino 
fluxes; the cross section for the ^He(''He,7)^Be reaction is the most important uncertainty for the 
calculated ^Be neutrino flux. 



Thispaper is part of a series that spans more than 40 
years The goals of this series are to provide increas- 
ingly more precise theoretical calculations of the solar 
neutrino fluxes and detection rates and to make increas- 
ingly more comprehensive evaluations of the uncertain- 
ties in the predictions. We describe here two steps for- 
ward (improved accuracy of the equation of state of the 
solar interior and some of the nuclear fusion data) and 
one step backward (increased systematic uncertainties in 
the determination of the surface composition of the Sun) . 

Using recent improvements in input data, we calculate 
the best-estimates, and especially the uncertainties, in 
the solar model predictions of solar neutrino fluxes. We 
compare the calculated neutrino fluxes with their mea- 
sured values. We stress the need for more accurate mea- 
surements of the surface composition of the Sim and of 
specific nuclear reaction rates. 

TableQlpresents, in the second (third) column, labeled 
BP04 (BP04-I-), our best solar model calculations for the 
neutrino fluxes. The uncertainties are given in column 2. 
BP04-I- was calculated with new input data for the equa- 
tion of state 0, nuclear physics [1,01, and solar composi- 
tion 0. BP04, our currently preferred model, is the same 
as BP04-I- except that BP04 does not include the most 
recent analyses of the solar surface composition Q , which 
conflict with helioseismological measurements. The error 
estimates, which are the same for BP04, BP04-I-, and ^^N 
(see Table nj, include the recent composition analyses. 

For the BP04 solar model, the base (mass) of the con- 
vective zone is O.TlSi?© (0.024 Mq), the surface heavy 
element to hydrogen ratio by mass, Z/X = 0.0229, the 
surface helium abundance is 0.243, and 1.6% of the lu- 
minosity is from CNO reactions. The central temper- 
ature, helium abundance, and Z/X, are, respectively, 
15.72 X 10*3 K, 0.640, and 0.0583. All of these values 
are in the acceptable range as determined by helioseis- 
mology. However, for BP04-f , the base of the convective 
zone (CZ) is Rcz/Rq — 0.726, which conflicts with the 
measured value of 0.713±0.001 (or, ±0.003, see Ref. 0). 
By examining a series of models, we have determined that 
the reason for the too-shallow CZ in the BP04-f model 



TABLE I: Predicted solar neutrino fluxes from solar mod- 
els. The table presents the predicted fluxes, in units of 
10^°(pp), lO^(^Be), 10*(pep,"N,i5Q)^ 10«(»B,"F), and 
lO^ihep) cm"^s"\ Columns 2-4 show BP04, BP04-I-, and 
our previous best model BPOO Columns 5-7 present the 
calculated fluxes for solar models that differ from BPOO by 
an improvement in one set of input data: nuclear fusion cross 
sections (column 5), equation of state for the solar interior 
(column 6), and surface chemical composition for the Sun 
(column 7). Column 8 uses the same input data as for BP04 
except for a recent report of the -I- p fusion cross sec- 
tion. References to the improved input data are given in 
the text. We use OPAL radiative opacities calculated for 
each chemical composition. The last two rows ignore neu- 
trino oscillations and present for the chlorine and gallium so- 
lar neutrino experiments the capture rates in SNU (1 SNU 
equals 10""^^ events per target atom per sec). Due to oscilla- 
tions, the measured rates are smaller: 2.6 ± 0.2 and 69 ± 4, 
respectively. We use the neutrino absorption cross sections 
and their uncertainties that are given in Ref. 



Source BPQ4 BP04-H BPOO Nucl EOS Comp 



pp 


5.94(1 ±0.01) 


5.99 


5.95 


5.94 


5.95 


6.00 


5.98 


pep 


1.40(1 ±0.02) 


1.42 


1.40 


1.40 


1.40 


1.42 


1.42 


hep 


7.88(1 ±0.16) 


8.04 


9.24 


7.88 


9.23 


9.44 


7.93 


'Be 


4.86(1 ±0.12) 


4.65 


4.77 


4.84 


4.79 


4.56 


4.86 


«B 


5.79(1 ±0.23) 


5.26 


5.05 


5.77 


5.08 


4.62 


5.74 




5.71(1 


4.06 


5.48 


5.69 


5.51 


3.88 


3.23 


15q 


5.03(1 t«:«) 


3.54 


4.80 


5.01 


4.82 


3.36 


2.54 




5.91(1 1^1) 


3.97 


5.63 


5.88 


5.66 


3.77 


5.85 



CI 


8.5t{i 


7.7 


7.6 


8.5 


7.6 


6.9 


8.2 


Ga 


13111^ 


126 


128 


130 


129 


123 


127 



is the lower heavy element abundance, ZjX — 0.0176. 
Therefore, we prefer BP04. 

The measurements from different solar neutrino exper- 
iments and the KamLAND reactor data ^3 can be 
combined in a global analysis to obtain the best empiri- 
cal values for the p-p, *B, and ^Be solar neutrino fluxes. 
We use the fluxes from the global analysis of Ref. 
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which allows all the solar neutrino fluxes to be free pa- 
rameters subject only to the luminosity constraint (i.e., 
energy conservation). Comparing the measured values 
with the theoretical predictions, we find for BP04: 

measured 

= (1.02 ± 0.02 ± 0.01)(/.(pp)thoory (1) 
0(^B)nioasurcd = (0.88 ± 0.04 ± O.23)0(^B) 

theory (2) 

</>(^Be)„,easurcd = (0.9lt^ j^ ± O.ll)0(^Be)thcory (3) 

In Eq. I^-Eq. 0, the la experimental uncertainties are 
given before the la theoretical uncertainties. 

The measured and theoretical values for the fluxes 
agree within their combined la uncertainties. The mea- 
surement error of the *B neutrino flux is smaller than the 
uncertainty in the theoretical calculation, but the oppo- 
site is true for the p-p and ^Be neutrino fluxes. 

Column four of Table presents the fluxes calculated 
using our previous best solar model, BPOO 0. The BP04 
best-estimate neutrino fluxes and their uncertainties have 
not changed markedly from their BPOO values despite 
refinements in input parameters. The only exception is 
the CNO flux uncertainties which have almost doubled 
due to the larger systematic uncertainty in the surface 
chemical composition estimated in this paper. 

We describe improvements in the input data relative 
to BPOO. Quantities that are not discussed here are the 
same as for BPOO. Each class of improvement is repre- 
sented by a separate column, columns 5-7, in Table 

Column 5 contains the fluxes computed for a solar 
model that is identical to BPOO except that we have 
used improved values for direct measurements of the 
''Be(p,7)^B cross section, 5*20 kGv(''Be -I- p) — 20.6 ± 
0.8 eV b ^, and the calculated p-p, So{pp) = 3.94(1 ± 
0.004) X 10-25 MeV b, and hep, So (hep) = (8.6 ± 1.3) x 
10-2°keV b, cross sections |^. The reactions that pro- 
duce the *B and hep neutrinos are rare; changes in their 
production cross sections only affect, respectively, the ^B 
and hep fluxes. The 15% increase in the calculated ^B 
neutrino flux, which is primarily due to a more accu- 
rate cross section for ^Be(p,7)^B, is the only significant 
change in the best-estimate fluxes. 

The fluxes in Column 6 were calculated using a refined 
equation of state, which includes relativistic corrections 
and a more accurate treatment of molecules The 
equation of state improvements between 1996 and 2001, 
while significant in some regions of parameter space, 
change all the solar neutrino fluxes by less than 1%. Solar 
neutrino calculations are insensitive to the present level 
of uncertainties in the equation of state. 

The most important changes in the astronomical data 
since BPOO result from new analyses of the surface chem- 
ical composition of the Sun. The input chemical compo- 
sition affects the radiative opacity and hence the physical 
characteristics of the solar model, and to a lesser extent 
the nuclear reaction rates. New values for C,N,0,Ne, 



and Ar have been derived [7| using three-dimensional 
rather than one-dimensional atmospheric models, includ- 
ing hydrodynamical effects, and paying particular atten- 
tion to uncertainties in atomic data and observational 
spectra. The new abundance estimates, together with 
the previous best-estimates for other solar surface abun- 
dances imply Z/X = 0.0176, much less than the 
previous value of Z/X — 0.0229 jl^J. Column 7 gives the 
fluxes calculated for this new composition mixture. The 
largest change in the neutrino fluxes for the p-p chain is 
the 9% decrease in the predicted ®B neutrino flux. The 
N and O fluxes are decreased by much more, ~ 35%, 
because they reflect directly the inferred C and O abun- 
dances. 

The CNO nuclear reaction rates are less well deter- 
mined than the rates for the more important (in the 
Sun) p-p reactions 0]. The rate for ^^N(p,7)^50 is 
poorly known, but important for calculating CNO neu- 
trino fluxes. Extrapolating to the low energies relevant 
for solar fusion introduces a large uncertainty. Column 
8 gives the neutrino fluxes calculated with input data 
identical to BP04 except for the cross section factor 
S'o(i'*N-hp) = 1.77 ± 0.2 keV b that is about hah the 
current best-estimate; this value assumes a particular R- 
matrix flt to the experimental data |l4l| . The p-p cycle 
fluxes are changed by only 1%, but the ^"^N and ^^O 
neutrino fluxes are reduced by 40% — 50% relative to 
the BP04 predictions. CNO nuclear reactions contribute 
1.6% of the solar luminosity in the BP04 model and only 
0.8% in the model with a reduced S'o(^'*N -I- p). 

Table ITU shows the individual contributions to the flux 
uncertainties. These uncertainties flx the accuracy to 
which a given input parameter should be determined. 
Improvements in input data enable more precise tests of 
stellar evolution via solar neutrino experiments. 

Columns 2-5 present the fractional uncertainties from 
the nuclear reactions whose measurement errors are most 
important for calculating neutrino fluxes. Unless stated 
otherwise, we have used throughout this paper the uncer- 
tainties estimated in Ref. 13] for nuclear cross sections. 

The measured rate of the '^He-'^IIe reaction, which after 
the inception of this series 0] changed by a factor of 4, 
and the measured rate of the ^Be -I- p reaction, which 
for most of this series has been the dominant uncertainty 
in predicting the ^B neutrino flux, are by now very well 
determined. If the published systematic uncertainties for 
the ■^He-'^He and ^Be -I- p reactions are correct, then their 
uncertainties no longer contribute in a crucial way to 
the calculated theoretical uncertainties (see column 2 and 
column 4 of Table ITTjl. These important improvements in 
our knowledge of low energy nuclear cross section factors 
have been achieved by increasingly precise, but always 
enormously difficult and beautiful, experiments carried 
out over four decades by many researchers in laboratories 
throughout the world. 

The most important nuclear physics uncertainty in cal- 
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TABLE II: Principal sources of uncertainties in calculating 
solar neutrino fluxes. Columns 2-5 present the fractional un- 
certainties in the neutrino fluxes from laboratory measure- 
ments of, respectively, the "^He-^He, ^He-''He, p-'^Be, and p- 
^''N nuclear fusion reactions. The last four columns, 6-9, give, 
respectively, the fractional uncertainties due to the calculated 
radiative opacity, the calculated rate of element diffusion, the 
measured solar luminosity, and the measured heavy element 
to hydrogen ratio. 



Source 3-3 3-4 1-7 1-14 Opac Diff L© Z/X 

pp 0.002 0.005 0.000 0.002 0.003 0.003 0.003 0.010 

pep 0.003 0.007 0.000 0.002 0.005 0.004 0.003 0.020 

hep 0.024 0.007 0.000 0.001 0.011 0.007 0.000 0.026 

■'Be 0.023 0.080 0.000 0.000 0.028 0.018 0.014 0.080 

*B 0.021 0.075 0.038 0.001 0.052 0.040 0.028 0.200 

0.001 0.004 0.000 0.118 0.033 0.051 0.021 0.332 

0.001 0.004 0.000 0.143 0.041 0.055 0.024 0.375 

0.001 0.004 0.000 0.001 0.043 0.057 0.026 0.391 



culating solar neutrino fluxes is now the rate of the '^He- 
''He reaction (column 3 of Table ITl)l . The systematic un- 
certainty in the the rate of ^He(^He, 7)^Be reaction(see 
Ref. [13j) causes an 8% uncertainty in the prediction of 
both the ^Be and the solar neutrino fluxes. It is as- 
tonishing that there has not been any progress in the past 
15 years in measuring this rate more accurately. 

For ^'*N(p,7)^^0, we have continued to use in Table UTI 
the uncertainty given in Ref. [l3< |. although the recent 
reevaluation in Ref. 0| suggests that the uncertainty 
could be somewhat larger (see column 7 of Table QJ. 

The uncertainties due to the calculated radiative opac- 
ity and element diffusion, as well as the measured solar 
luminosity (columns 6-8 of Table lll| . are all moderate, 
non-negligible but not dominant. For the ^B and CNO 
neutrino fluxes, the uncertainties that are due to the ra- 
diative opacity, diffusion coefficient, and solar luminosity 
are all in the range 2% to 6%. 

The surface composition of the Sun is the most prob- 
lematic and important source of uncertainties. System- 
atic errors dominate: the effects of line blending, de- 
partures from local thermodynamic equilibrium, and de- 
tails of the model of the solar atmosphere. We assume 
that the uncertainty in all important abundances is ap- 
proximately the same. We have defined previously the 
3(7 range of Z/X as the spread over all modern de- 
terminations (see Refs. P, y, ^3)' which now implies 
A{Z/X)/{Z/X) = 0.15 (Icr), 2.5 times larger than the 
uncertainty adopted in BPOO. The most recent uncer- 
tainty quoted for oxygen, the most abundant heavy ele- 
ment in the Sun, is similar: 12% 7]. 

Heavier elements like Fe affect the radiative opacity 
and hence the neutrino fluxes more strongly than the 
relatively light elements 0. This is the reason why the 



difference between the fluxes calculated with BP04 and 
BP04-I- (or between BPOO and Comp, see TablcQJ is less 
than would be expected for the 23% decrease in Z/X. 
The abundances that have changed significantly since 
BPOO (C, N, O, Ne, Ar) are aU for fighter species for 
which meteoritic data are not available. 

The dominant uncertainty listed in Table ITU for the ^B 
and CNO neutrinos is the chemical composition, repre- 
sented by Z/X (see column 9). The uncertainty ranges 
from 20% for the **B neutrino flux to - 35% for the CNO 
neutrino fluxes. Since the publication of BPOO, the best 
published estimate for Z/X decreased by 4.3cr(BP00 un- 
certainty) and the estimated uncertainty due to Z/X in- 
creased for ^B (^^O) neutrinos by a factor of 2.5 (2.8). 
Over the past three decades, the changes have almost 
always been toward a smaller Z/X. The monotonic- 
ity is surprising since different sources of improvements 
have caused successive changes. Nevertheless, since the 
changes are monotonic, the uncertainty estimated from 
the historical record is large. 

We list below our principal conclusions and their impli- 
cations. First, the experimentally-determined values for 
the p-p, ''Be, and ^B solar neutrino fluxes are in agree- 
ment with the values predicted by the standard solar 
model. In each case, the agreement is accidentally better 
than the combined Icr uncertainties. More precise mea- 
surements of the p-p and ^Be fluxes will test critically the 
theory of energy generation in the solar interior. Second, 
recent precise measurements or improved calculations of 
nuclear reaction parameters, the equation of state, and 
the surface chemical composition of the Sun have re- 
fined the input data to the solar model calculations but 
have not changed the calculated neutrino fluxes outside 
the previously-quoted theoretical Icr uncertainties (see 
columns 3-6 of Table |I| and Ref. [3). Third, the rate of 
the reaction ^He(^He,7)^Be is the largest nuclear physics 
contributor to the uncertainties in the solar model predic- 
tions of the neutrino fluxes in the p-p chain. In the past 
15 years, no one has remeasured this rate; it should be 
the highest priority for nuclear astrophysicists. For the 
important ^Be neutrino flux that can be measured in the 
BOREXINO [ill and KamLAND ^ detectors, there is 
currently an 8% uncertainty due to the cross section for 
■^He(^He,7)^Be. Fourth, the cross section for the reaction 
^'*N(p,7)^^0 is the largest nuclear physics contributor to 
the uncertainties in the calculated CNO neutrino fluxes. 
It is important to measure this cross section more accu- 
rately in order to understand well energy production in 
stars heavier than the Sun. Neutrino oscillation stud- 
ies can use the fluxes from both the low ^''N model and 
BP04 (see TablePi; the results for oscillation parameters 
should be essentially identical. Fifth, the largest uncer- 
tainty in calculating all the solar neutrino fluxes is now 
the uncertainty in the measured surface composition of 
the Sun (see Table lll|) . Unfortimately, the principal un- 
certainties in inferring the composition are systematic. 
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The uncertainties could be estimated more reliably, and 
perhaps reduced, by comparing the results from a variety 
of different measurements of the same elemental abun- 
dances and by different groups of astrophysicists making 
comparably sophisticated two and three dimensional hy- 
drodynamical models of the solar outer region. Sixth, 
the recent reanalyses of the solar chemical composi- 
tion imply a lower surface heavy element abundance and 
consequently a base of the convective zone that conflicts 
with helioseismological measurements 0. For this rea- 
son, we have not used the most recent low heavy element 
abundances in our preferred model, BP04. The low heavy 
element abundances lead in BP04-I- to a slightly better 
agreement with the neutrino measurement, but the 
improvement between best-estimates is smaller than the 
la uncertainty. 

We make three recommendations for future work, 
based on the known dependences of neutrino fluxes on 
input parameters First, the low energy cross sec- 

tion of the ■^He(^He,7)^Be reaction should be measured 
to better than ±5% {la){a factor of two improvement) 
in order that the uncertainty in this reaction not limit 
the interpretation of future ^Be solar neutrino experi- 
ments. Second, the uncertainty in the surface heavy el- 
ement abundances (particularly elements like iron that 
contribute most significantly to the radiative opacity 0) 
should be reduced to less than ±0.02 dex (a factor of 
3 improvement) in order that the calculational uncer- 
tainty from the composition not exceed the current error 
(±4% 0, 0, 0|) in the empirical determination of the 
^B neutrino flux. Third, the uncertainty in the low en- 
ergy extrapolation of the rate of the -'^"'N(p,7)^'''0 reaction 
must be < 25% in order that the p-p flux can be used 
for a precision measurement of the neutrino mixing angle 
012 [Uj and an accurate test of stellar evolution theory. 
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